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Abstract: As a main component in water balance, evapotranspiration (ET) is of great importance for 
water saving, especially in arid and semi-arid areas. In this study, the FAO (Food and Agriculture 
Organization) Penman-Monteith model was used to estimate the magnitude and temporal dynamics of 
reference evapotranspiration (ETo) in 2014 in subalpine meadows of the Qilian Mountains, Northwest 
China. Meanwhile, actual ET (ETà was also investigated by the eddy covariance (EC) system. Results 
indicated that ET. estimated by the EC System was 583 mm, lower than ETo (923 mm) estimated by the 
PAO Penman-Monteith model in 2014. Moreover, ETo began to increase in March and reached the peak 
value in August and then declined in September, however, ET. began to increase from April and reached 
the peak value in July, and then declined in August. Total ET. and ET» values during the growing season 
(from May to September) were 441 and 666 mm, respectively, which accounted for 75.73% of annual 
cumulative ET, and 72.34% of annual cumulative ET, respectively. A crop coefficient (&.) was also 
estimated for calculating the ET,, and average value of & during the growing season was 0.81 (ranging 
from 0.45 to 1.16). Air temperature (Ta), wind speed (m), net radiation (R,) and soil temperature (T) at the 
depth of 5 cm and aboveground biomass were critical factors for affecting &, furthermore, a daily 
empirical & equation including these main driving factors was developed. Our result demonstrated that 
the ET, value estimated by the data of & and ETo was validated and consistent with the growing season 
data in 2015 and 2016. 
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1 Introduction 


Water balance strongly influences the structure, function, and productivity of terrestrial 
ecosystems (Burman and Pochop, 1994). A critical component of water balance is actual 
evapotranspiration (ET,) (Mitchell et al., 2009; Valipour et al., 2017), and determination of water 
balance depends on precise estimation of ET; (Jensen et al., 1990; Allen et al., 1994; Jensen et al., 
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1997; Allen et al., 1998; ASCE-EWRI, 2005; Veroustraete et al., 2008; Allen et al., 2011; Semalty 
et al., 2011; Hssaine et al., 2018). Understanding of ET, is also important for water utilization and 
management, especially in arid and semi-arid regions (Li et al., 2014; Zhao and Zhao, 2014). 
Qilian Mountains, located in Northwest China, are water sources for many inland rivers in Hexi 
Corridor. Subalpine meadows (grasslands) are a major land use in these mountains where heavy 
grazing not only impacts grass growth conditions but also affects ET. Therefore, it is imperative 
to estimate ET, of these alpine meadows with consideration of grazing impact on ETÀ. 

A number of methods and models have been developed and used to estimate ET, such as 
Bowen (Ohmura, 1982; Sun et al., 2006; Fang et al., 2018), Lysimeter (Howell et al., 1991; Gao 
et al., 2015), eddy covariance (EC) (Swinbank, 1951; Massman, 2000), water balance (Mastrorilli 
et al., 1998; Saadi et al., 2018) and FAO (Food and Agriculture Organization) Penman-Monteith 
(P-M) model (Allen et al., 1998). This study focused on two methods: EC and FAO 
Penman-Monteith model. The EC is a micro-meteorological approach widely used to study ET, at 
the ecosystem level (Massman, 2000; Falge et al., 2001; Aubinet et al., 2012). The FAO 
Penman-Monteith model is a modified version of P-M model (Monteith, 1965) recommended by 
the FAO to estimate ET, for various land uses including grassland (Allen et al., 1998). For each 
land use, ET, is estimated using reference evapotranspiration (ETo) values which in turn depend 
on land use specific crop coefficient (ke) (Allen et al., 1998; Yang et al., 2013; Khoshravesh et al., 
2017). The k. values, for a specific land use, are generally assigned by the growth stage of the 
plants species. For example, the ko values of 0.40, 1.05, and 0.85 in grasslands are respectively 
reported for initial-, middle-, and late-season growth stages of grass (Allen et al., 1998). Studies 
also show that some additional factors, such as meteorological data and vegetation biomass, also 
affect ke value (Lockwood, 1999; Williams and Ayars, 2005; Zhou and Zhou, 2009; Yang and 
Zhou, 2011; Zhang et al., 2012). Therefore, accurately calculating ke can provide more reliable 
estimation of ET.. 

The ET, for subalpine meadows in the Qilian Mountains has been estimated (ranging from 2.88 
to 3.26 mm/d) using micro-lysimeter method (Song et al., 2004; Yang et al., 2013). However, the 
ET, estimation by micro-lysimeter is influenced by the weather conditions, i.e., it might be 
over-estimated during sunny days and under-estimated during rainy weather (Zhang et al., 2012; 
Zhou and Zhou, 2009). Thus, it is unsuitable for long-term ET. estimation during all weather 
conditions. To fill this knowledge gap for grasslands in the Qilian Mountains, we used (1) EC 
system to observe ET, at daily time scale, and (2) FAO Penman-Monteith model to estimate ETo. 
The objectives of this study are (1) to estimate ET. based on k: and ETo values of the subalpine 
grasslands in the Qilian Mountains, and (2) to develop and validate a model of k. based on 
meteorological factors and vegetation biomass of subalpine grasslands in the Qilian Mountains. 


2 Materials and methods 


2.1 Study area 


The study area is located in Tianlaochi Catchment in the upper reaches of the Heihe River 
(38°23'56"—38°26'47"N, 99°53'57"—99°57'10"E) with the elevation ranging from 2660 to 4400 m 
a.s.l. (Fig. 1). Experimental site for this study is located at 3070 m altitude in subalpine meadow 
in the Qilian Mountains. The study area belongs to a typical alpine semi-arid climate, 
characterized by long and cold winters, and short and warm summers. The mean annual 
precipitation is about 435.0 (+5.6) mm with 84.2% precipitation occurring during May to 
September. The annual mean temperature is about —0.6°C—2.0°C, and mean January and mean 
July temperatures are —13.1°C and 12.1°C, respectively. The soil type is mountain grey cinnamon 
with up to 50 cm depth and gravel below 50 cm. Vegetation community mostly belongs to annual 
grasses, i.e., Elymusnutans griseb and Foeniculum vulgare, and their roots grow up to 30 cm. 


2.2 Data collection 


2.2.1 Meteorological and soil physical data 


At the study site, a weather station and an EC system were set up in 2011 and 2013, respectively. 
Table 1 shows the list of meteorological factors for which data are being collected at half an hour 
intervals and soil physical data at the study site. 
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Fig. 1 Geographic location of the study area and distribution of observation instruments in the experimental site. 
EC, eddy covariance. 


Table 1 Summary ofthe meteorological and soil physical data observed at the weather station 


Parameter Company Country Observational elevation Abbreviation 
Net radiation Gene USA 2m Ry 

Short wave radiation Gene USA 2m R; 

Soil heat flux Gene USA —5 cm G 

Wind speed Hobo USA 2m 

Air temperature Hobo USA 2m Ta 

Relative humidity Hobo USA 2m RH 
Pressure Hobo USA 2m P; 
Precipitation Hobo USA 1.8m P 

Soil temperature Decagon USA 5, -10, —20, and —50 cm Tss, Ts10, Ts20, and Tiso 
Soil water content Decagon USA 5, -10, —20, and -50 cm SWCs, SWCjo, SWC20, and SWCso 


2.2.2 Flux measurements 


The EC system was used to measure CO: fluxes, sensible heat (H), and latent heat (LE) at 2.2 m 
above the ground surface from 2013. A CO2/H2O infrared analyzer (Li-7550; LI-COR, Inc., 
Lincoln, NE, USA) and a three-dimensional supersonic anemometer (CSAT-3; Campbell 
Scientific, Inc., UT, USA) were mounted on a horizontal bar extending from a tower. The 
CO>/H20 sensor was installed in down-wind direction of supersonic anemometers. The CO2/H20 
analyzer was calibrated after every six months. Short wave radiation (Rs) from the sky and long 
wave radiation (Rı) from the land surface were measured with a net radiometer (CNR-1; Kipp and 
Zonen, Delft, The Netherlands) at 1.5 m above the ground surface. The net radiation (Rn) was 
calculated from the difference of Rs and Rj. Soil heat fluxes (G) were recorded by a net radiometer 
at 10 cm soil depth. The Ry was partitioned into H, LE, and G. 
The energy balance ratio (r) was calculated using the following equation (Gu et al., 1999): 


r=(ġ (LE + H))/ ÈR, -G)). (1) 


The item LE+H measured by the EC system seemed to be underestimated since the average 
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value of r was about 0.75 for the entire observation period from 2014 to 2016, which fell in the 
median region of reported energy closures (0.55—0.99) (Wilson et al., 2002). The lack of energy 
balance closure had also been reported by Stannard et al., (1994) and Aubinet et al. (2012), and 
energy balance closure had become accepted as an important new test of EC (Anderson et al., 
1984; Mahrt, 1998). We were not trying to specify a particular cause for the imbalance because 
several possibilities may be involved in the lack of energy closure (Wilson et al., 2002). 

In long-term field measurements, gaps in data are a known issue. The gaps in data occur either 
due to non-recording (instrument malfunction/power failure) or rejection of certain data that do 
not meet quality control requirements. The quality of flux data collected for this study was 
assessed through a widely adopted flag system (Mauder and Foken, 2004), which lead to discard 
some data recorded at night time. The data gap filling was therefore, performed using a lookup 
table approach that considered both covariation of fluxes with meteorological variables and the 
temporal autocorrelation of the fluxes (Falge et al., 2001; Reichstein et al., 2005). 

2.2.3 Biotic data 

At the experimental site, six quadrats (0.5 mx0.5 m) were randomly selected to collect 
aboveground biomass. Biomass was measured every three days by cutting the plants at ground 
level and loading the plants into sampling bags during the growing season (May—October). Each 
sample was numbered and dried at 65°C until a constant weight was obtained. 


2.3 Methods 


2.3.1 FAO Penman-Monteith model 


The FAO Penman-Monteith model (Allen et al., 1998) was used to calculate ETo, which is 
expressed as follows: 


ee (R,-G)+— 20 Jun. (2) 
A\A+Y y +A\T+275 

ET, =k, ET, ’ (3 ) 

y =y(1+0.33u), (4) 


where À is the water latent heat of vaporization (MJ/kg); A is the slope of the vapor pressure 
versus temperature curve (kPa/°C); R, is the net radiation (MJ/(m?°-d)); G is the soil heat flux 
(MJ/(m*.d)); y is the psychrometric constant (kPa/°C); T is the daily mean temperature at 2 m 
above the ground surface (°C); u is the wind speed at 2 m above the ground surface (m/s); D is the 
vapor pressure deficit (kPa); ETcm is the estimated actual evapotranspiration (mm/d) using the 
FAO-P-M model; and ks is the crop coefficient. 

2.3.2 EC model 

EC model was used to calculate the ET.: 

ET.=1000LE/1 pw, (5) 
where LE is the latent heat flux (MJ/(m?.d)); A is the water latent heat of vaporization (MJ/kg); 
and pẹ is the water density (kg/m’). 

2.3.3 Crop coefficient (ke) 

The daily kc was obtained by dividing ET, values with ETo values for the year 2014. The obtained 
k: was used to develop two models: one with meteorological factors only, and the other using both 
meteorological factors and biotic factors. Both models were developed using the data of 2014, 
and validated using the data of 2015 and 2016. 

2.3.4 Evaluation methods 

In the study, the linear correlation coefficient (R°), slope, Relative Root-Mean-Squared Error 
(RRMSE), Coefficient of Determination (CD) and Index of Agreement (IA) were used to examine 
the extent of the differences between the calculated and observed values of ET.. The slope was 
the relative between the straight line and horizontal axis. The R* value was used to assess the 
correlation between the prediction results with actual situation (Gao et al., 2013). The RRMSE 
was a measure of relative magnitude of residuals; and the models with smaller RRMSE value are 
desirable (Gao et al., 2015). The IA was a measure to evaluate the correlation between the 
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observed and simulated values, and IA values closer to 1 are desirable (Li et al., 2014). The CD 
was used to measure dispersion of simulated values from the mean observed values, and CD 
greater than 0.8 indicates that the result was acceptable (Cai et al., 2007). The R? and slope were 
calculated in Sigmaplot 11.0 (Alexandris and Kerkides, 2003; Yang and Zhou, 2011). The 
mathematical expressions of RRMSE, CD, and IA were presented as follows: 


n E 2 
nse = 2a E), (6) 
O 


210-9) | D 


n 2 
"i dial - 9) (8) 
2a 
where O; is the observed value; O is the mean observed value; E; is the simulated value; and n is 
the number of samples. 
2.3.5 Vapor pressure deficit (VPD) 


VPD is the difference between the pressure of water vapor held in saturated air at a given 
temperature and the pressure of water vapor that exists in the ambient air. It is calculated as 
follows: 


17.27xT 


VPD = 0.61008 x e7 #2373 x (1 — 5, (9) 


where T is the daily mean temperature at 2 m above the ground surface (°C); RH is the relative 
humidity. 


3 Results 


3.1 Variations of meteorological and soil physical factors 


Temporal variations in meteorological factors are shown in Figures 2 and 3. The mean daily 
values of Rn, Rs, Ta, VPD, and u were 12.43 MJ/(m?.d), 22.40 MJ/(m?.d), —0.57°C, 0.658 kPa and 
1.21 m/s, respectively. During the growing season, the sequence of SWC was 
SWC20>SWCs>SWC10>SWCso and that of T; (soil temperature) was Tss>Ts10>Ts20>Tss0. Both Ta 
and Rn reached the minimum in December and the maximum in July. VPD reached the maximum 
of about 1.180 kPa during the growing season, whereas u reached the maximum in February and 
the minimum in July. 


3.2 Monthly variations of ETo and ET. 


Figure 4a shows that annual ETo was 923 mm, which was higher than the annual ET, (583 mm) in 
2014. ETo gradually increased from January and reached the peak value in August and then 
declined in September. Whereas, ET, began to greatly increase in April and reached the peak level 
in July and then it started to decline in August. The total ET, and ET» values during the growing 
season (from May to September) were 441 and 666 mm, respectively, which accounted for 
75.73% and 72.34% of annual ET, and ETo, respectively. 

Figure 4b compares the temporal variation of ET. and precipitation (P). The annual values of 
ET, and P were 583 and 546 mm, respectively. The annual ratio of ET. to P is 1.067. However, 
the values of ET. and P during the growing season (from May to September) were 441 and 540 
mm, respectively, and the ratio was 0.818. 


3.3 Temporal variations of ke 


As shown in Figure 5, ke quickly increased from May to June, and remained a high level in July 
and August. It declined in September and was about 0.2 at the end of October. During the growing 
season, kc exhibited a mean value of 0.81 and fluctuated within the ranges of 0.45-1.16. 
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Fig. 2 Temporal variations of meteorological and soil physical factors from 2014 to 2016. (a) net radiation (Rn) 
and short wave radiation (Rs); (b) soil temperature (Ts) at the depths of 5, 10, 20, and 50 cm (i.e., Tss, Ts10, Ts20 and 
Ts50, respectively); (c) soil water content (SWC) at the depths of 5, 10, 20, and 50 cm (i.e., SWCs, SWC10, SWC20 
and SWC50, respectively). 
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Fig. 4 Annual variations in actual evapotranspiration (ET-) and reference evapotranspiration (ETo). (a) annual 
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Fig. 5 Variation of crop coefficient (ko) at 10-d intervals in 2014. Bars indicate standard errors. 
3.4 Factors impacting ke 


3.4.1 Effect of environmental factors on ke 


In the study, the Pearson's product-moment correlation coefficients between k: and environmental 
factors were computed at a daily scale (Table 2). The significant correlations (P<0.01) of ke with u, 
jm Ta, Rn and Tss suggest that these factors are critical for controlling ke. The relationships of ke with 
— P, SWCs, Ts10 and 7,20 were also significant (P<0.05). 

l The significant negative correlation between u and RH (Table 2) suggests the role of wind in 
the blowing away water vapours. The SWCs mainly come from precipitation, i.e., the infiltration 
into soil layer through rainfall, so there are good correlations between precipitation and SWC. 
Solar radiation provides energy to soil that in turn increases soil temperature. Using significant 
climatic factors (Table 2), we developed the following regression model to estimate daily ke: 

k, = 0.650 + 0.224u — 0.0387, —0.030R, — 0.002 P + 0.570SWC, 


(10) 
+0.1107.; — 0.0167., —0.0097.,,. 
According to Equations 3 and 9, we calculated the daily ETemoai as follows: 
ET. oa = K, X ET, = (0.650 + 0.224u — 0.0387, —0.030R, — 0.002 P + 0.570SWC, (11) 


+0.1107,, -0.0167,,, 0.0097.) xETy, 


where ETemoai is the ET. calculated by Equation 11. 
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Table 2 Pearson's product-moment correlation coefficients between daily crop coefficient (ko) and 
meteorological and soil physical factors, i.e., wind speed (u), air temperature (Ta), relative humidity (RH), net 
radiation (Rn), air pressure (Pa), precipitation (P), soil water content (SWC) at the depths of 5, 10, 20 and 50 cm, 
and soil temperature (7s) at the depths of 5, 10, 20 and 50 cm during the growing season in 2014 


ke u Ta RH Ra P, P SWCs SWCio SWC2  SWCso Tss Ts10 Ts20 Tss0 
ke 1.000 
u 0.194" 1.000 
Ta 0.177 0.099 1.000 
RH —0.109  -0.899"* —0.019 1.000 
Rn 0.512 0.037 0.203" -0.012 1.000 
P, —0.139* 0.043 —0.090 -0.165" —0.039 1.000 
P -0.055  -0.356"" 0.108 0.394** 0.169** -0.169"* 1.000 
SWCs —0.049* -0.542*" -0.045 0.600 -0.120° -0.141™ 0.155%" 1.000 


SWC10 0.028  -0.307" 0.077 0.383** -0.004 —0.257** -0.044 0.751™* 1.000 
SWC 0.084 -0.079 0.072 0.151* -0.005 -0.153™" -0.047 0.565°" 0.741" 1.000 
SWCs0 0.122 -0.185" 0.470°* 0.292** 0.107 -0.166"" 0.106 0.582°* 0.662" 0.666%" 1.000 


Tss 0.173 —0.069 0.933" 0.189** 0.259** —0.197*™* 0.237** 0.082 0.164" 0.148" 0.595** 1.000 

Ts10 0.150% -0.184 0.891™* 0.306"" 0.237 -0.201™ 0.280" 0.192" 0.241" 0.218" 0.679"* 0.980" 1.000 

Ts20 0.119* —0.320*™ 0.794** 0.442% 0.180"" —0.171*™* 0.265™* 0.365* 0.375" 0.296** 0.778"* 0.901°* 0.962" 1.000 

Ts50 0.091  -0.403"" 0.651" 0.514" 0.067 -0.119"*  0.141** 0.530" 0.514" 0.342™* 0.792™* 0.744" 0.831" 0.942" 1.000 


Note: ** and * indicate significance at 0.99 and 0.95 confidence levels, respectively. 


3.4.2 Effects of biotic factors on ke 


Biomass directly reflects the conditions of a biological community, so we selected it as a biotic 
factor. The results showed that the biomass gradually increased in May, rapidly increased in June 
and July, and reached a peak in August (with a maximum of 505 g/m’), and then began to 
decrease rapidly in September (Fig. 6). 
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Fig.6 Variation of biomass during the growing season in 2014. Bars indicate standard errors. 


The regression analysis showed a significant positive correlation between biomass and ke (Fig. 
7). The linear regression equation is as follows: 
k,=0.0011biomass+0.5687. (12) 
Considering both meteorological and biotic factors, we calculated daily ETcmoa2 as follows: 
ET cmoa2=(0.650+0.224u—0.038 7,—0.030R,—0.002P+0.570S WCs5+0.1107ss—0.0167510— 
0.009 7520) X (0.0011 X biomass+0.5687)xETo, (13) 


where ET cmoa2 is the ET, calculated by Equation 13. 
3.5 Testing estimated EC, 


We validated Equations 11 (ETemodi) and 13 (ETcmoaz2) using the EC system data collected during 
the growing seasons in 2015 and 2016. The estimated values were close to the observed values (Fig. 
8). The slope, R? and IA of ETemoa2 were greater than those of ETcmoai (Table 3), implying that the 
correlation between observed values and simulated values by ETemoa2 is better than that by 


201906.00047v1 


chinaXiv 


ChinaX iv (ERAT 


GAO Yunfei et al.: Actual evapotranspiration of subalpine meadows in the Qilian Mountains... 


y=0.001 1x+0.5687 
2= (1.6900 
jals P<0.01 


0 50 100 150 200 250 300 350 400 450 500 550 
Biomass (g/m’) 


Fig. 7 Scatter plot between crop coefficient (kc) and biomass 


ETemoai. The RRMSE obtained by ETcmoa2 was smaller than that by ETcmoai, suggesting a higher 
accuracy of ETemoa2 than that of ETcemoai. The values of CD for both ETemoai and ETemoa2 were 
greater than 0.8, indicating both regression equations are good. Overall, the performance of 


ETcmoa2 that used both meteorological and biotic factors is better than that of ETcmoa: (only 
meteorological factors). 
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Fig. 8 Temporal variations of measured ET (ETc) and simulated ET (ETemoa1). (a) ETc and ETecmoai in 2015; (b) 
ET: and ETemoa2 in 2015; (c) ETc: and ETemoa! in 2016; and (d) ET, and ETemoa2 in 2016. 
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Table 3 Summary of the parameters to assessing the accuracy of the two models 


Slope R RRMSE IA CD 
Variable 
2015 2016 2015 2016 2015 2016 2015 2016 2015 2016 
ET emodi 0.875 0.835 0.709 0.867 0.077 0.102 0.930 0.874 1.002 0.812 
ET omoa2 0.921 0.890 0.866 0.906 0.050 0.092 0.972 0.916 0.942 0.835 


Note: R?, correlation coefficient; RRMSE, Relative Root-Mean-Squared Error; IA, Index of Agreement; CD, Coefficient of 
Determination. 


4 Discussion 


4.1 Variations of ET», ET, and P 


The result shows that ET. was lower than ETo. The maximum of ETo appeared in August, and the 
minimum was in December. However, the maximum of ET. appeared in July and the minimum 
occurred in December. 

According to Equation 2, the meteorological and soil factors were the only factors that affected 
ETo. However, the ET, was affected not only by meteorological factors but also by biotic factors 
such as biomass, when Rs, Rn, Ta, and RH increased from April and reached the peak values in 
August and then declined (Figs. 2 and 3). ET. began to increase in April and reached the 
maximum in July with the canopy development during the growing season. 

The annual ratio of ET, to P (i.e., 1.067) indicates that P was less than ET, (Hou et al., 2011). 
However, the ratio of ET, to P during the growing season (from May to September) was 0.818. 
This suggests that adequate water availability supports the plant growth in this period. The value 
of ETo during the growing season was 666 mm, and the ratio of ETo to P during the growing 
season was 1.240. This result was similar to the researches in some grassland ecosystems (Wever 
et al., 2002; Li et al., 2005). However, our result was higher than a non-irrigated pasture (Sumner 
and Jacobs, 2005) and an alpine winter pasture (Li et al., 2014). 

ET process is governed by energy exchange at the vegetation surface and is limited by the 
amount of energy available (Semalty et al., 2011). ETo in our study area is caused by the climatic 
conditions. Rn in the study area was higher than that in the Qinghai-Tibetan Plateau grassland (Gu 
et al., 2008), so ETo was higher than that reported by Li et al. (2014) who conducted their study in 
the Qinghai-Tibetan Plateau grassland. Apart from energy enhancing the ETo in some extent, the 
high SWC and P could offer water for the need of plants (Yang and Zhou, 2011; Simona and Rita, 
2013; Li et al., 2014). High 7; may inspire plant activity and make plant absorb more water (Gao 
et al., 2015), while low RH may improve the water vapor exchange between the atmosphere and 
vegetation (Gao et al., 2015). 


4.2 Temporal variations of ke 


ko quickly increased from May to June, remained at a high level in July and August, and then 
declined in September. The range of ke in the study (0.45—1.16) was larger than that of ke reported 
in FAO 56 (0.30-1.05) (Allen et al., 1998), and it was higher than that of k: observed in an alpine 
winter pasture (0.30—0.92) (Li et al., 2014), in a typical steppe (0.32—0.68) (Miao et al., 2009) and 
in a temperate desert steppe of Inner Mongolia (0.02-0.50) (Yang and Zhou, 2011; Zhang et al., 
2012). Our result indirectly showed that the subalpine meadows in the Qilian Mountains featured 
more adequate moisture and better plant growth conditions than the other semi-arid and arid 
grassland ecosystems. 


4.3 Factors impacting k, 


The highest correlation between kc and environmental factors suggested that Rn was the most 
important factor for determining ke. Similar findings have been reported for an alpine winter 
pasture (Li et al., 2014). The higher Rn made a significant contribution to the ke (Gao et al., 2015). 
However, Yang and Zhou (2011) selected a temperate desert steppe as their study area, which 
water was a key environmental factor pointed out that SWC was the most important factor, and 
Hou (2010) stated temperature was the most important meteorological factor. These differences 
may result from the different environmental conditions of plant growth. 
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In addition to meteorological factors, biotic factors can greatly affect ke (Tyagi et al., 2000; 
Williams et al., 2003; De Medeiros et al., 2005; Li et al., 2014). Biomass directly reflects the 
conditions of a biological community, so we selected it as a biotic factor to estimate ke. 

Figure 9 shows that the EC calculated by both Equations 11 and 13 is underestimated in 2015 
and 2016. FAO Penman-Monteith model performs well under the abundant water condition. The 
underestimation of ETo and ET., in our study area is likely due to the erratic precipitation which 
is not sufficient to meet water requirement criteria of FAO Penman-Monteith model. 
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Fig. 9 Scatter plots between ETc and ETcmoai in 2015 (a), ETc and ETemoa2 in 2015 (b), ETc and ETcmoai in 2016 
(c), and ET. and ETcmoa2 in 2016 (d) 


4.4 Testing estimated EC, 


Equations 11 (ETemoai) and 13 (ETcmoa2) were validated using the EC system data collected during 
the growing season in 2015. The simulated accuracy by the Equations 11 and 13 was compared 
using slope, R?, RRMSE, IA, and CD. 

The values of RRMSE were smaller, and the results of the model correlation were better. In 
Table 3, the slope, R? and IA of ETcmoai were greater and the RRMSE was smaller than those of 
ETcmoa, and the CD values of ETemoai and ETcmoa2 were greater than 0.8. The results indicate that 
the fitting method that comprehensively considers meteorological factors and vegetation 
conditions can be used to derive simulated ET, results that were closer to the actual observed ET, 
results. 


5 Conclusions 


The FAO Penman-Monteith model was used to simulate the ETo over the subalpine meadows in 
the upper reaches of the Heihe River. Meanwhile, the ETc was measured by EC system. During 
the growing season, the ratio of ET; to ETo is 0.663, and the ratio of ET; to P is 0.818. Both ratios 
explained the habit of meadows in the study area is good. The k. was used to estimate the ETc. 
The average value of k: during the growing season was 0.81 (ranging from 0.45 to 1.16). Rn, Ta, u, 
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and 7;5 were the critical factors for controlling ke. Pa, SWCs, Tsio and T3290 also showed a 
significant relation to ke. Finally, we developed a daily empirical k. equation based on biomass, Ta, u, 
Rn, P, SWCs, Tss, Ts10, and Ts20 to estimate the daily ke, and used the k: and ETo to estimate the 
ET.. The model validation with 2015 and 2016 growing season data shows that the daily 
empirical k: equation performed well in this study area. 
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